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SUMMARY 



The initial surface reactions of the extrinsic coagulation pathway on live cell membranes were 
examined under flow conditions. Generation of fXa (activated coagulation factor X) was measured 
on spherical monolayers of epithelial cells with a total surface area of 41-47 cm 2 expressing TF (tissue 
factor) at > 25 fmol/cm 2 . Concentrations of reactants and product were monitored as a function of 
time with radiolabeled proteins and a chromogenic substrate at resolutions of 2-8 s. At physiological 
concentrations of fVIIa and fX, the reaction rate was 3.05 ± 0.75 fmol fXa/s/cm 2 , independent 
of flux, and 10 times slower than that expected for collision-limited reactions. Rates were also 
independent of surface fVIIa concentrations within the range 0.6-25 fmol/cm 2 . The transit time of 
fX activated on the reaction chamber was prolonged relative to transit times of nonreacting tracers 
or preformed fXa. Membrane reactions were modeled using a set of nonlinear kinetic equations and 
a lagged normal density curve to track the expected surface concentration of reactants for various 
hypothetical reaction mechanisms. The experimental results were theoretically predicted only when 
the models used a slow intermediate reaction step, consistent with surface diffusion. These results 
provide evidence that the transfer of substrate within the membrane is rate-limiting in the kinetic 
mechanisms leading to initiation of blood coagulation by the tissue factor pathway. 



INTRODUCTION 



Blood coagulation reactions mediate fibrin deposition in hemostasis and many pathological pro- 
cesses. Blood clots are directly implicated in the lethal complications of cardiovascular disease and 
contribute significantly to the pathogenesis of infectious, autoimmune, and neoplastic diseases (1-6). 

The blood coagulation process is initiated by an assembly of complexes comprised of an essential 
cofactor, TF (tissue factor) and a protease component, fVIIa. The functional complex, TF/fVIIa, 
cleaves the natural substrates, fVII, fIX (factor IX), and fX at specific sites, generating, fVIIa, fIXa, 
and fXa, respectively (4-6). Factors VII, IX, and X circulate in the blood and extravascular fluids 
(7-10), while TF is expressed on the membranes of many extravascular tissues (|ll|)- The anatomic 
distribution of cells expressing TF is consistent with its role as the initiator of hemostatic reactions. 
Cell surfaces in contact with blood do not appear to express functional TF constitutively. However, 
inflammatory stimuli induce expression of functional TF on endothelial cell membranes and blood 
monocytes (12-14). 

Factors VII, IX, and X are vitamin K-dependent proteins, and their functional interaction with 
negatively charged procoagulant membranes has a calcium-dependent, electrostatic component (15- 
19). The interaction sites are located in highly homologous 7— carboxyglutamic acid (Gla) rich 
regions near the N-terminus of all vitamin K-dependent coagulation proteins (4, 19). The specific 
binding and functional kinetics of interaction between coagulation proteins and biological mem- 
branes have been studied extensively under equilibrium steady-state conditions (20-27). While 
equilibrium binding parameters vary significantly among vitamin K-dependent proteins, adsorp- 
tion parameters are similar, suggesting nonspecific initial contact (28, 29). Anionic phospholipid 
membranes modify the apparent kinetic parameters of coagulation reactions relative to kinetics in 
solution. The membrane effect is manifested by a large decrease in the apparent K m of substrates to 
values far below their respective plasma concentrations. The mechanisms by which this effect man- 
ifests itself during TF-mediated coagulation remain speculative. Achieving useful time resolutions 
has been one of the main obstacles to developing experimental systems to study the presteady-state 
transients of coagulation factor adsorption and activation on cell membranes. Blood clotting in vivo 
and in vitro can be completed faster than the sampling intervals of traditional batch systems used 
to measure membrane reactants. 

Measurements of fVIIa binding and fXa generation on intact cell membranes under steady-state 
conditions indicate that TF/fVIIa functional activity is fully expressed before the binding interaction 
between fVIIa and TF reaches equilibrium (24, 25). Furthermore, under steady-state conditions, 
the overall rate of coagulation substrate activation on membranes pre-equilibrated with enzyme was 
close to the theoretical collisional limit (|27|) . These findings suggest intermediate noncovalent steps 
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on the membrane linking the initial adsorption step to the assembly and catalysis of substrate in 
the A-E-S (activator-enzyme-substrate). 

Coagulation zymogens and active proteases are subject to local microcirculation controls (3, 
7-10, 14) since they are found in extravascular lymphatic, synovial, and alveolar fluids. The im- 
portance of flow control in coagulation reactions has been demonstrated in vivo. Tracer studies 
with radiolabeled fibrinogen and vasoactive agents indicate a direct correlation between changes in 
vascular permeability and fibrin deposition (|3|). Several studies using lipid-coated capillaries also 
indicate that flow rates influence the activity of coagulation proteases (30, 31). In the present study, 
we use high-resolution tracer-dilution analyses (32-35) along with numerical modeling to identify 
the surface and flow-dependent kinetics of fX activation via the TF pathway. We show that the 
generation of fXa from plasma fX proceeds via an intermediate step within the membrane. For 
reactions initiated with fVIIa and fX, the rate of this step and of the overall reaction is limited by 
the transfer of fX from the adsorption sites to the catalytic sites on the cell's surface. 

EXPERIMENTAL PROCEDURES 

Cell Culture and Reaction Chambers-Cells and cell cultures have been described and character- 
ized in detail elsewhere (29, 34). Briefly, Vero cells (American Tissue Type Collection) were grown 
to confluency on microcarrier beads (Cytodex 2, Pharmacia) of 150 fim average diameter. Reaction 
chambers were assembled with a 2-3 ml suspension of cell-covered microspheres loosely packed in 
a thermoregulated column fitted with flow adapters. A schematic of the reaction flow chamber is 
shown in Figure |[ 

Cell viability and metabolic integrity, as demonstrated by amino acid uptake, were maintained 
for periods exceeding the kinetic measurements described here (0). The cells expressed surface- 
TF const it utively and interacted functionally with human fVIIa and fX with apparent steady-state 
kinetic parameters similar to other procoagulant cells, as reported previously (29, 35). The TF 
activity of the monolayer suspension in the reaction chamber was equivalent to 6149 ± 847 fmol/ml 
of recombinant TF reconstituted into phospholipid vesicles. The chambers were perfused with 
medium (M-199, Gibco), buffered with HEPES, and supplemented with 0.1 mM of ovalbumin and 
3 mM CaCl2- The overall geometrical parameters and TF content in the reaction chamber are 
summarized in Table I. 

Tracer- dilution techniques-Concentrations of reactants and products in the bulk, flowing, aque- 
ous phase and on the cell surface during fXa-generating reactions were measured using double and 
triple-tracer techniques. These techniques were adapted from previously described and validated 
methods in perfused organs (32-34, 36, 37). Coagulation initiation reactions were carried out at 
physiological concentrations, under flowing conditions. Reactions were initiated under flow by in- 
jecting 48-500 ng of fVIIa and 3000-9000 ng of fX in a 1.23 ml reaction chamber []. The injected 
reactants are quickly dispersed in the reaction chamber and attain maximal initial concentrations 
of approximately 0.8-8 nM fVIIa and 45-135 nM fX. concentrations but diminish over time. Control 
tracer 14 C-labeled ovalbumin (833 ng) was also included in the bolus injection yielding an initial 
maximal concentration of approximately 14 nM. To measure membrane concentrations, 3 H-labeled 
fVIIa was used an adsorbing tracer. Reactants and control tracer were in a perfusing medium 
containing 0.1 mM unlabelled ovalbumin and 3 mM CaC£2- 

After introducing reactants into the chamber via the inflow tubing, the effluent was collected 
in 72 samples of 46 ± 2.8 /d each at time resolutions ranging from 2-10 s per sample, depending on 
the perfusion rate. At flow rates of 13 /d/s, the typical flow velocity in the reaction chamber was 
estimated at < 0.2 mm/s. This value is within the ranges measured in rabbit microvasculature in 
vivo (|38|). 

Samples were collected in microtiter wells preloaded with TRIS buffer solution, pH 8.3, con- 
taining 0.2 M EDTA and 0.8 M NaCl. A final 5 ml sample was collected to complete the recovery 
of tracer and to determine perfusion rates. Standard curves for 14 C and 3 H tracers were con- 
structed from serial dilutions of the solutions used in the bolus injection. Concentrations of control 

1 These amounts of reactant were rapidly (1-2 s) injected in bolus of 100 jii 
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and test radioactive tracers in standard dilutions and effluent samples were measured by scintil- 
lation counting. The concentration of product, fXa, was measured by amidolytic assays (13, 14) 
with chromogenic substrate (Methoxycarbonyl-D-cyclohexylglycil-arginine-p-n itroanilide-acetate) 
before scintillation counting of radioactive tracers. 

Functional tests for TF Activity in Ce//s-Functional activity of TF was determined from fXa 
generation rates in purified systems using amidolytic assay and recombinant proteins as standard. 
Recombinant TF used to construct standard curves was relipidated into 30:70 PS/PC (phosphaty- 
dylserine/phosphatydyl choline) vesicles as before (17, 31). The molar ratio of TF/lipid in the 
standard TF-PS/PC preparations was 1/4100. The TF activity measured in intact monolayers was 
68.5 ± 19% of that measured in monolayers lysed by freezing/thawing. 

Radioactive Tracers-The control tracer used to measure concentrations in bulk aqueous phase 
was 14 C-labeled ovalbumin (Sigma, St. Louis, MO) with a specific activity of 33 /zCi/mg. The test 
tracer for adsorption measurements was fVIIa radiolabeled with tritium using the technique of Van 



Lenten and Ashwell (f40|) , with modifications (fill). Labeled preparations had specific activities of 



1.8x10 s cpm/mg of fVIIa and a functional activity comparable to unlabeled factor in clotting tests 
and activating mixtures with purified components (p9|). 

Reaction scheme and mathematical model-The surface reactions leading to fXa generation were 
analyzed according to the following scheme: 

k+E 

fVIIa + TF £ E 
E + fX ^ a E-fX^fXa + E, 

k — n 



where E = fVIIa-TF is the fVIIa and TF complex ("enzyme") that forms and dissociates with rate 
constants k + E,k-E, respectively. The substrate-enzyme complex denoted by E-fX associates and 
dissociates with a second-order rate-constant, k +a , and a first-order rate-constant, k_ a , respectively. 
The effective rate of product (fXa) formation from the complex and its irreversible release are 
denoted by the first-order rate-constant, k + . Denoting the surface concentrations of each species 
by (r 1; T 2 , r 3 , r 4 , T 5 , T 6 ) = (fVIIa, fX, TF, E, E-fX,fXa), the full kinetic equations consistent with 
flD are 

f i = -k +E TiT 3 + k_ E ^i - (3^! + ctiCi(t) 
f 2 = -k +a T 2 T 4 + k_ a T 5 - (3 2 T 2 + a 2 C 2 (t) 
T 3 = -k +E ^i^z + k_ E ^4 

(2) 

r 4 = k+ET^-3 - fe_ B r 4 - A; +a r2r 4 + (k_ a + k + )r 5 



f 5 = A; +a r 2 r 4 - (k^ a + k + )T 5 

f 6 = fc+r 5 + a 6 c 6 - /3 6 r 6 , 

where Ti(t) = dTi(t)/dt, with i = (1, 2, 3, 4, 5, 6). The constants k±E and k± a correspond to effective 
rates of TF and fX binding interactions with TF and E respectively. These coefficients include the 
time delays of all intermediary processes on the membranes before the interactions. The time 
distributions of these unspecified processes are accounted for in our numerical predictions of the 
time course for the overall fX activating process. The possibility of inhibition or fX/fXa-destroying 
sinks is precluded from our data since, within experimental error, all the absorbed fX is recovered 
as fXa. 
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In the above nonlinear differential equations, /3, and aij are desorption and adsorption rates, 
respectively. Since under our experimental conditions the total area fraction of adsorbed species is 
negligible, species adsorption from bulk is simply proportional to the bulk concentration, Cj(t), at 
the surface of each microsphere. As the fluid passes through the ensemble of microcarriers, certain 
flow lines are faster or slower than the mean flow velocity, resulting in a distribution of reactant 
velocities. A "lagged normal density curve" (LNDC) has been successfully used to approximate 
the dispersion resulting from the combined effects of random velocity distribution and molecular 
diffusion in the human circulatory system (36, 37). We find good agreement between a fitted 
lagged density curve and the sequentially measured concentrations in the outflow of the reaction 
chamber (Fig. |3|). Therefore, to simplify the modeling process, we assume that the dispersion and 
diffusion of all species are equal and use the LNDC to approximate the source, Ci(t), surrounding 
each microcarrier. The parameters used in the lagged density curve will reflect chamber packing 
characteristics, bulk diffusion constants, and the imposed volume flow rate, Jy. Additional details, 
analysis, and simplifications of Eqs. (0) are provided in the Appendix. 

Calculation of reactant concentrations in membranes-The proportion of fX and fVIIa adsorbed 
from the flowing phase into the membrane was determined from the difference between the nor- 
malized concentrations of control, 14 C-labeled ovalbumin, and 3 H-labeled fVIIa. Concentrations of 
factor Vila adsorbed at time t were estimated using 

Ut) * {[ u C(t)] - [ 3 R(t)]) Q t St\ (3) 

where [ 14 C] and [ 3 H] are the fraction of the total nonadsorbed control and adsorbed test tracer, 
respectively, collected in the effluent. Qr{in fmol) is the total amount of fVIIa added, and St is 
the total membrane surface area (~ 41-47 cm 2 ). Previous studies measuring adsorption of various 
coagulation factors, including fVIIa and fX, indicated that adsorption rates are proportional to their 
aqueous-phase concentration and not significantly different among vitamin K-dependent proteins 
(28, 29). Based on these data, the normalized concentration of 3 H-labeled fVIIa was used to trace 
both fVIIa and fX adsorption. For substrate fX, the membrane concentration r 2 (t) was determined 
by 

r a (f) « ([ 14 C(t)j - [ 3 H(i)] - [fXa(t)]) Q T Sj\ (4) 

where Qt is the total amount (fmol) of fX added, and [fXa(t)] is the fraction of that total released 
into the aqueous phase as fXa. 

Miscellaneous-Tissue factor antigen expressed by the cells was determined in cell lysates using 
a commercial ELISA kit with recombinant soluble TF as standard (American Diagnostics). Pro- 
tein determinations were performed in the same cell lysates with a commercial reagent (Biorad 
Laboratories), using bovine serum albumin as standard. Coagulation fVIIa was human recombi- 
nant, kindly donated by Dr. Ulla Hedner (Novo Nordisk, Denmark). Recombinant tissue factor 
(used as standard in ELISA and functional TF determination in cells) was purchased from Amer- 
ican Diagnostics. Human fX and fXa were purchased from Enzyme Laboratories. Data reduction, 
plotting, and statistical analyses were performed using Stat View software (Brain Power, Inc.). Nu- 
merical solution of the nonlinear kinetic equations (0) was performed using adaptive Runge-Kutta 
implemented through Matlab. 

RESULTS 

Rate of fXa generation at different fVIIa concentrations-The generation of fXa from fX on live 
procoagulant cell membranes was examined in reaction chambers filled with spherical cell monolay- 
ers. The geometrical and flow characteristics of these reaction chambers are summarized in Figure 
|I] and Table I, respectively. The distribution of concentrations of reactants in the flowing bulk 
aqueous phase was followed using control tracer 14 C ovalbumin. Reactions were initiated with fX 
and fVIIa, and the product, fXa, was measured by amido lytic assay. Reactions were followed until 
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70-90% of the nonreacting control tracer was recovered in the effluent. The amounts of 14 C and 3 H 
tracer collected and the amount of 3 H adsorbed to the cell are shown as functions of time in Figure 
@A. 

The time evolution of reactants inside the reaction chamber can be fairly well approximated 
by the LNDC, as shown in Figure || This agreement indicates that dispersion of reactants in 
the chamber due to the random flow distribution and diffusion is qualitatively similar to that 
encountered in human circulation (36-38, 43, 46-47). Figure [3] also illustrates the time/concentration 
distribution of fXa generated in the 1.23 ml reaction chamber in a typical reaction initiated with 
500 ng of fVIIa and 9000 ng of fX. 

Under these conditions, aqueous-phase concentrations ranged from 0.3-10 ±2 nM and from 4- 
137 ± 29 nM for fVIIa and fX, respectively. A time trace measuring the total amount of fXa collected 
and fX adsorbed on cell membranes is shown in Figure ^|B. Factor Xa profiles were weakly sigmoidal 
with a linear middle segment. Table II shows that the average production rate, calculated from the 
linear segment, did not change when average concentrations of fVIIa in aqueous phase were decreased 
by ten-fold, from 5 nM to 0.5 nM (membrane concentrations ranged from 0.7-25.0 fmol/cm 2 ). No 
fXa was generated in the absence of fVIIa and reaction rates did not differ significantly when fVII 
was substituted for fVIIa. The observation that maximal constant catalytic activity is reached at 
very low concentrations of TF-fVIIa complexes allows for simplifying substitutions in the model 
equations (Eqs. fj) for Ti and T 4 (Appendix). 

Comparison between the reaction rate and the theoretical collisional rate-The independence of 
average reaction rates and enzyme concentrations suggests that substrate transfer to the catalytic 
sites is rate limiting. Two possibilities were investigated: (i) the rate-limiting step may occur during 
the adsorption of reactants from bulk to the cell membranes, and (ii) the rate-limiting step occurs 
after the adsorption step. To differentiate these two possibilities, the rate of fXa generation was 
compared to the theoretical collisional rate between reactants and microspheres, given the aqueous 
phase concentrations of fX used and the flow rates, Jy, imposed. Since average adsorption rates 
were shown to approach or to exceed the collisional limit (^), activation rates below this limit 
support the second possibility. 

Theoretical steady-state collisional rates were calculated from the aqueous-phase concentrations 
of fX and the radius of the spherical microcarriers using Smoluchowski's relationship for steady-state 
diffusion (42, 57) 

k co u « D h2 C h2 (t)R-\ (5) 

where k co u is the collision rate between reactant molecules and a unit area of membrane (collisions/cm 2 /s), 
Di^ is the diffusion constant for fVIIa, fX in water (~ 5 x 10~ 7 cm 2 /s), R ~ 7.5 x 10~ 3 cm is the 
microcarrier radius, and C\p is the fVIIa, fX concentration (molecules/cm 3 ). 

Figure ^ contrasts the number of fXa molecules released by the monolayer and the aqueous 
phase concentration, C 2 (t), of fX as a function of time. Since the collisional rate follows Eq. 0, 
collision-limited rates are expected to be directly proportional to C^{€). However, the rate of fX 
activation on the monolayer was not correlated with fX-membrane collisions. The rate of fXa 
production (molecules/cm 2 /s) reached maximal values after the peak in C^i) and collisional rates. 
Furthermore, high fXa rates were sustained during the rapid decrease in collisions between fX and 
the membrane, following the concentration peak. Averaged over 13 experiments, the activation 
rate was 3.05 ± 0.72 fmol/cm 2 /s, corresponding to 1.8 ± 0.43 xlO 9 (molecules/cm 2 /s), below the 
theoretical maximum of 2.4 ± 0.57x 10 10 (collisions/cm 2 /s). 

Apparent second-order rate coefficients, calculated from initial reaction rates and aqueous-phase 
concentrations of fVIIa and fX, did not have a constant value but increased continuously during 
the observation time. Using aqueous-phase fX and membrane-phase fVIIa to calculate apparent 
second order rate coefficients also resulted in increasing coefficient values, consistent with the obser- 
vation that the average reaction rates are essentially independent of fVIIa concentration (Table II). 
These results indicate that the initial adsorption from bulk to membrane is not rate-limiting in the 
overall fXa generation reaction. The results also imply that fX is activated via a membrane-bound 
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intermediate rather than directly from the bulk aqueous phase. 

Mean transit times of fX through the reaction chamber-The existence of a slow membrane step 
was further investigated by comparing average transit times of the fXa generated in the chamber to 
the transit times of preformed fXa. The presence of a rate-limiting step between membrane adsorp- 
tion and catalytic cleavage is expected to delay the transit of the substrate that is adsorbed and 
catalyzed as compared to bulk aqueous-phase reactants. A mean transit time, To, was determined 
from the concentration of fXa and control tracer in 72 consecutive samples of the effluent according 
to the expression 

where [*C(t n )] is the fraction (of the total amount added or theoretical maximum) of either control 
tracer or fXa collected in aliquot n at time t n ~ n x (2 — 8) sec, depending on the particular 
experiment. Results shown in Table III indicate that the mean transit time of fX, To(fX), activated 
in the reaction chamber is increased relative to the To of aqueous-phase control tracer. In contrast, 
To(fXa) for fXa formed before being introduced in the reaction chamber is indistinguishable from 
that of the control tracer. Furthermore, the increase in To(fX) is inversely correlated with flow 
rate. These results are consistent with a slow membrane step following the fast, flow-dependent 
adsorption step. 

Kinetic modeling of surface reactions-The hypothesis that the reaction pathway proceeds with 
fast equilibration of enzyme activity followed by a rate-limiting step involving reactant surface dif- 
fusion was also tested by comparing experimental measurements to the solutions of the kinetic 
equations (0). Equations (0) were solved numerically using initial estimates for intrinsic rate con- 
stants based upon results of previous steady-state kinetic studies (15, 20, 21, 25-27). Heuristic 
arguments for initial guesses for all the rate parameters are provided in the Appendix. A continu- 
ous function for aqueous-phase concentrations Ci(t) is derived from a least-squares fit to a LNDC 
(Appendix), shown in Fig. |3|. The remaining parameters in the model were then adjusted until 
the best visual fit of r 6 (t) to fXa collected was achieved. Since T 2 (t) was indirectly measured and 
subject to larger experimental errors, we only varied rate parameters to get an order-of-magnitude 
agreement between the measured fX (Fig. ^) an d T 2 (t) (Fig. f|B), using measured fXa and r 6 (t) to 
more precisely fit the parameters. The solutions and the associated best-fit parameters are shown 
in Fig. |J We found that the magnitudes of Ti(t) match the measurements only when the amount of 
TF assumed in the simulations was 0.32 fmol/cm 2 , much smaller than the actual amount expressed 
on the cell membranes. This finding is consistent with our hypothesis that enzyme complexes form 
domains, further developed in the Discussion Section. 

The model also shows that within reasonable ranges, the shape and magnitude of the product 
curve, r 6 (t), are sensitive to k +a ,k + and aij, but less sensitive to the other parameters. If the 
association step, k +a , were fast, the theoretical model would predict a premature overproduction of 
fXa, as shown in Figure |5|A. The sensitivity to a slow intermediate membrane step associated with 
k +a is shown in Fig. |5|C, while the corresponding predicted values for effective enzyme, TF, and fX 
on the membranes are shown in Fig. [5p. Note also that for parameters differing from those used 
in Fig. ||, the magnitudes of r 2 (t) change dramatically and are no longer close to the measured fX 
(Fig.|). 

For the parameters used to fit the measurements in Fig. 2, the numerical solution for E (r.^) 
plateaus to a value ~ T\ after t ~ 40 s and remains nearly constant for the duration of the 300 
s interval under consideration. This quasi-steady-state exists even for the cases where k +a is too 
large (Figure |5[B) or too small (Figure |5p). We show in the Appendix that this quasi-steady-state 
behavior allows us to define an approximate effective rate constant (s _1 ) 

k + k +a T* A 

k '* = k^TkZ (7) 

that approximately describes the rate of fXa production on the cell membranes via 
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T 6 Rj k e ffY 2 {t) + adsorption/desorption terms. 



(8) 



For large fc + (the fast chemical step), k +a becomes the limiting rate, since k e ff ~ k +a . The 
experimental data are consistent with model predictions both qualitatively and quantitatively, when 
k +a is in the range expected for lateral diffusion of proteins on membranes. An estimate for a 
mechanistically relevant diffusion length can be derived from 



where D% is the surface diffusion constant of fX in the cell membranes. Using a typical value 
-Qsurf ^ ^Q-io cm 2^ s (pjj^ anc } fc^g ^ 0.02 s _1 obtained from the mathematical model, we find that 
Id ~ 0.7/xm. However, can be shorter, if obstructions in the membrane hinder surface diffusion 
and reduce D s ^ ur . 

Factor Xa generation rate as a function of surface density of reactants-The observations de- 
scribed in the previous sections indicate that reaction rates are not directly related to the aqueous- 
phase concentration of substrate. To further investigate the rate-limiting step we analyzed reaction 
rates as a function of flow rates and membrane concentrations of substrate. 

The instantaneous fraction of membrane fX converted to fXa was not directly proportional to 
the fX concentration on the membrane. Instead, at all flow rates tested, it increased linearly with 
time. Interestingly, the increase in the proportion of adsorbed fX encountering catalytic sites per 
unit time was essentially independent of flow. This observation suggests that adsorbed fX is not 
immediately available to the catalytic sites and is consistent with a rate-limiting surface diffusion 
process. That the fraction of membrane fX available for catalysis increases with time suggests 
diffusive transfer of fX from initial adsorption sites to the catalytic sites. The value of the slope of 
the function T 6 (t)/T 2 (t) was 0.001 - 0.0005 and largely independent of flow rate (Table IV). 

Table IV also lists the product yield and average reaction rates measured at different flow rates. 
The yield was strongly correlated with flow rate (correlation coefficient 0.86), while the average 
reaction rate was independent of flow rate. Again, these results are as expected for a kinetic 
mechanism that includes a flow-dependent adsorption step followed by a rate-limiting intermediate 
step on the membrane. 



In this paper, we have analyzed the surface reaction kinetics of blood coagulation initiated by 
interaction among coagulation factors on biological membranes. Reactions were initiated on live 
epithelial cells expressing TF with physiological concentrations of fVII/fVIIa and fX. The aqueous 
and membrane concentration of reactants as well as product formation were measured at time 
resolutions relevant to plasma clotting. Whereas adsorption of vitamin K-dependent proteins on 
procoagulant membrane surfaces has been shown to be fast and correlated to aqueous-phase flux 
(2"9D, the rate of fXa generation was independent of both enzyme density on the membrane and 
flow rate, Jy. Moreover, using tracer dilution analyses we found that the transit time of the fX 
participating in the reaction was prolonged relative to transit times of nonreacting control tracers. 
Rate coefficients calculated from reaction rates and either the aqueous or membrane concentration 
of reactants changed with time. Flow velocities influenced the total amount of reactants adsorbed 
to the membrane and the total yield, but not the intrinsic rate of product formation. Taken 
together, the experimental results provide evidence for a kinetic mechanism limited by a slow 
transfer of substrate between initial membrane adsorption sites and reaction sites. The hypothesis 
of a slow membrane step was further tested by numerically solving a set of nonlinear kinetic equations 
describing the evolution of all membrane reaction species. The experimental results were reproduced 
when the rate-limiting step followed the adsorption of substrate to the membrane and preceded the 




(9) 



DISCUSSION 



8 



chemical catalysis. Testing the alternative hypothesis of either slow adsorption or slow catalytic 
steps resulted in product yield and profiles markedly different from those observed experimentally. 

Maximal surface catalytic activity was observed when the membrane concentration of fVIIa was 
at least two orders of magnitude lower than the surface concentration of TF, estimated either by 
immunoassay or functional tests. The experimentally measured fXa levels best matched solutions 
of the kinetic equations (Eqs. 0) when the intrinsic rate constants shown in Fig. ^ were used along 
with a maximal enzyme concentration of T4 ~ 0.32 fmol/cm 2 . This observation suggests that only 
a fraction of the available membrane E = fVIIa-TF is involved in catalysis as would be expected if 
catalytic sites were in large molar excess over the substrate. However, the substrate concentration 
reached, > 100 fmol/cm 2 , is much higher than TF or enzyme concentrations. 

These results can be explained using a model in which the enzyme concentration, T^, is in relative 
local excess, and only a fraction of the surface enzyme effectively participates in catalysis. If enzyme 
complexes fVIIa-TF form domains, the fast catalytic cleavage reaction would be expected to occur 
only near the perimeter of these enzyme domains, where the substrate initially encounters the 
enzyme after diffusing a typical distance, £d, following adsorption. The interiors of these domains 
are rarely accessed by the rapidly converted substrate and therefore do not participate in the overall 
reaction. This process is shown schematically in Figure || 

Surface segregation of molecules is a common phenomenon. It has been shown that even small 
molecules can form domains in lipid monolayers (44, 48, 49) and bilayer vesicles (0). A large 
body of experimental research also provides evidence that protein and lipid domains exist on live 
cell membranes (50, 51). For example, adhesion molecules localize at sites of cell-cell contact, 
and receptors are often found concentrated at the tips of filopodia and lamelopodia of moving 
cells. Neurotransmitter receptors in postsynaptic terminals have also been shown to form dynamic 
aggregates (53, 54). Direct visualization reveals localization of certain proteins in areas of high or 
low curvature in artificial vesicles ([55]). We have evidence of the existence of domain formation on the 
membrane of the epithelial cell line used in these studies. Using gold immunochemistry on cells fixed 
after short a exposure to fVIIa, the enzyme was localized primarily on the ruffled border of the cell 
membrane. Furthermore, analysis of nearest-neighbor distances indicated a nonrandom distribution 
of the enzyme (|29|). Although the mechanisms of domain formation are unknown, possibilities 
may involve electrostatic or dipole- induced phase transitions ([4*BD, and membrane elasticity-induced 
protein-protein attractions floTj|). 

The flow characteristics and time resolution achieved with this experimental system is relevant to 
reactions on biological membranes after exposure of TF to flowing plasma coagulation proteins. The 
experimental approaches and mathematical model used here to identify the early kinetic mechanisms 
of fXa generation will also be useful for studying novel pharmacokinetic mechanisms occuring on 
biomembranes. 

APPENDIX 

In this section, we give details of the mathematical model (Eqs. 0) and the associated approx- 
imations used in its analysis. Since no measurable amount of fXa is generated by fVIIa and fX in 
bulk solution, we have assumed that chemical reactions can only occur when molecules are adsorbed 
on the surface of each cell-covered sphere. 

The kinetic equations (Eqs. ||) are solved numerically using finite difference approximations. 
All surface densities, Tj, are in units of fmol/cm 2 , while all bulk concentrations, Cj, are measured 
in units of pmol/cm 3 . With this convention, the rates take on the following units: [j3] = [k- a ] = 
[k-E\ = [k+] = s _1 , [k + E\ = [k+ a ] = cm 2 /(fmol s), and [a] = 10~ 3 cm/s. 

We assume that the adsorption rates of species from the bulk onto the microsphere surfaces are 
proportional to the local bulk concentration. For the sake of completeness, and to motivate more 
quantitative modeling, we write the governing equations for surface adsorption of reactants. In the 
bulk phase, the concentration of species % follows the convection-diffusion equation: 

d t C + V ■ VC = DiV 2 C r>R, (10) 
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where C = C±,C2,Cq, and Di are the bulk solution concentrations (number per volume) and 
associated diffusion constants of fVIIa, fX, and fXa, respectively. Although a closure relation is 
required to specify the detailed velocity field, V, around each sphere, we will assume that the 
identical microspheres each feel an equivalent, averaged, effective flow velocity, V. The boundary 
conditions within continuum theory at the sphere surface are found by balancing the diffusive flux 
with the desorption and adsorption rates on the surface of each microsphere (at r = R), 



Did r Ci(r = R,t)= a^Ci(R,t) - frF^t), 



where the area fraction available for adsorption $ ~ 1 at low coverage. 

The above equations constitute the exact continuum equations for the species in solution. The 
complete set of equations consists of the convection-diffusion equation (Eq. [T0|) , the boundary 



2 



can be realized by assuming that 



conditions (Eqs. |TT|), and the surface reaction equations (Eqs. 

Significant simplifications and decoupling of some of Eqs. 
the bulk concentrations at the spheres' surfaces Ci(R,t) can be approximated by the LNDC fl37j). 
As the fluid passes through the ensemble of cell-covered microspheres, certain flow lines are faster 
or slower than the mean flow velocity. The reactants in the aqueous phase are randomly advected 
through the microsphere chamber at a distribution of velocities. The LNDC is a convolution of 
random advection velocities with molecular diffusion and has been used to approximate advection- 
diffusion in blood flow Q3"7|). Although the source concentration, Ci(t), may depend on the position 
of the microsphere within the reaction chamber, we assume that the dispersion and diffusion of all 
species are equal and use the lagged density curves to approximate the source, Ci(t), surrounding 
each microsphere. The parameters used in the LNDC depend upon average microsphere packing, 
the bulk diffusion constants, and the imposed constant volume flow rate, Jy. 



The concentration of the i species in a random flow environment is assumed to obey 



Tr 



d d(R,t) 
~~dt 



+ Ci(R,t) 



1000m,- 



exp 



T i )/a l 



(12) 



Jy(2^f 12 

where m« is the total number of femtomoles of species % added via the bolus injection, and Jy is 
the constant flow rate measured in fil/s. The intrinsic delay time, T«, is inversely related to the 
mean |V|, while a, and Tj describe the width and effects of molecular dispersion, respectively. The 
amount of spreading embodied in <Tj is proportional to the bulk diffusion constant, Di. The solution 
to the initial value problem (Eq. [12]) is 



Q(t) 



1000m 



27iJi 



V 



Erf 



Erf I 



Or 



(t ~ T- t 



\v2Tj V2cn 



(13) 



The concentration, Cj, above is given in nM units. These solutions determine the sources, 
ctiCi(R,t), for the surface kinetic equations (0). We take the entire reaction chamber and the inlet 
and outlet tubes to constitute a single, effective flow system. The zero used in Eq. [L3| corresponds 
to the time when nonbinding species are first detected (for the experiment in Fig. |3|, approximately 
36 s after adding reactants). Upon fitting (by adjusting <jj, Tj, Tj until a local minimum in the least- 
squares is found) the parameters in Eq. [13] to the concentration, we find (for this experiment at 
J v ~ 13.4 (A/s) that a* w 21 s, r » 91 s, and T » 19 s for % = 1, 2. The fitted LNDC is shown in 
Fig. |. 

The initial rate parameters used to solve Eqs. ^| were estimated as follows. The absorption 
rates of fVIIa and fX, from previous studies, were found to be similar and close to the collisional 
limit (|29|). Here, we first assumed that these adsorption rates were diffusion-limited. Such high 
absorption rates were needed to obtain the right magnitudes of fXa formation, regardless of the 
other rate parameters. Enough reactant must simply reach the membranes within the time limit 
imposed by the flow rate, Jy. The maximum rate of particles reaching and absorbing into the 
membrane of an isolated microsphere's surface, in the absence of flow, is given by 
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Jj < 4irRDiCi(r = oo, t). (14) 

This upper limit assumes that every molecule coming into contact with the sphere is absorbed. 
With a diffusion constant of Di ~ 7.5 x 10~ 7 cm 2 /s, the absorption rate under zero flow conditions 
is approximately «j < 1000Di/R ^0.1 cm/s (the factor 1000 converts pmol/cm 3 to fmol/cm 3 ). 
Now consider the effects of advection due to the imposed volume of flow, Jy. Purcell (|57|) gives an 
expression for the flux to a spherical surface under flow: 



RV 



1/3 



Ji(V) < toRDidit) ^—j . (15) 

Although at first glance the V 1 ^ 3 dependence is weak, quantitative changes in a due to flow can 
dramatically influence the yield of the surface reactions. For the density of microspheres and flow 
rates used, we estimated the typical velocity in the reaction chamber to be V ~ Jy/A e g « 0.02 cm/s, 
where Jy ~ 13.4/il/s, and the reaction chamber effective cross-sectional area A e g < 7r(0.7cm) 2 due 
to partial obstruction. Therefore, we used a\ 2 ~ 0-6 cm/s as an initial guess for the absorption 
rates in the reaction scheme (Eqs. §). 

For the association rates cti, we used the surface diffusion of fVIIa and fX to set upper limits. 
For a surface diffusion constant of ~ 10~ 10 cm/s, we found that k +a ~ k + E < 0.1cm 2 /(fmol 
s). Moreover, from previous aqueous-phase equilibrium binding studies, k_E/k+E ~ 10~ 10 M (24- 
26). To estimate the corresponding ratio for two-dimensional reactions, we made a qualitative 
estimate by assuming the energetics of E formation are not significantly different from those in 
bulk. Thus, 10~ 10 M corresponds to a typical particle-particle distance of 2.5 /im. Translating this 
to a surface density, we estimated very roughly that for the surface enzyme formation reaction, 
k-E/k + E ~ 10~ 2 fmol/cm 2 . Finally, from previous steady-state measurements on cell membranes 
at saturated concentrations of fX, k + > 14/s (pop. Within these limits, we explored the parameter 
space to obtain a reasonable fit to the data. The data and the fit of r 6 (t) calculated from (H) are 
shown in Fig. £|A. The concentration of product collected, Q^it), is calculated from 

Q e (t) « (3 6 T 6 (t)S T /J v , (16) 

where — 47 cm 2 is the total membrane area in the reaction chamber. The membrane concentra- 
tion r 2 is shown in Fig. |]B. 

The results were consistent with the data in Table II and showed that after a short initial 
transient, the concentration of membrane enzyme (r 4 ) reaches a plateau. This behavior permitted 
simplification of Eqs. |2| and approximate analytic solutions for the surface concentrations r 2 6 (t). 
For large k + + k_ a , the concentration T 5 (t) (E-fX) is always small. From a typical simulation Figs. 
[I]. [5] we observed a short transient in T 3 (TF). At times beyond this transient, T 3 w 0, and the 
enzyme concentration, T 4 , reaches a nearly steady value: 

r* « ^r 1 r 3 . (17) 

k-E 

From the equation for r 1; we see that over long times, fVIIa approximately follows the adsorption 
and desorption processes, 

fi»aiC7i(t)-Ari. (18) 

It is evident from Fig. f|-[5]that T 4 (E) also reaches a quasi-steady state shortly after TF. Therefore, 
setting r 4 w 0, 

1 - /, ' + " J : -r 2 = k eff v 2 . (19) 



k + + fc_ a 
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The remaining time- dependent surface quantities at these quasi-steady-state times obey 



+ 



a 2 C 2 (t) 
a 6 C 6 (t) 



(20) 



r 2 (t) \ = ( -keff o \f r 2 (f) 

T 6 (t) J V k eff -fa ) V r eW 

where /c ej y given by Eq. [5] is the effective rate of conversion from fX to fXa during quasi-steady-state 
times when the enzyme concentration is T 4 « r|. Assuming that OqCq is negligible, Eqs. ^0] admit 
analytic solutions and, considering the approximate nature of our model, a further simplification 
can be made: Using Eq. [16|, the second equation in (|20|) becomes 



Qe = k e fff3 6 



Si 



k 



eff 



-Qt 



(21) 



Therefore, an independent measurement of T 2 and the collected product Q 6 can be used to 
estimate the unknowns (3$ and k e g. Although in our analyses we have numerically solved the full 
kinetic equations (Eqs. 2), a simplified set of equations (Eqs. [H| and |2"0|) provide an analytic model 
to the reaction kinetics for times beyond the initial short transient (T> 40 s for the run shown in 
Fig. §. 
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Figure 1: Schematic of the reaction chamber and cell-covered microspheres. The reactive 
surface in the reactor is the surface of viable Vero cells grown to confluency on micro carriers with 
the indicated dimensions. The microcarriers (~ 5 — 20 x 10 4 ) are packed in a thermoregulated 
column fitted with flow adapters and perfused at constant flow rates of 5 — 25 /xl/s. Reactants and 
control tracer are added via the inflow (lower) and collected via the outflow (upper) tubing in 72-140 
consecutive samples of 46 ± 8.9 /xl each. For most experiments, reactants are added as a rapid bolus 
and reactions followed for 150-300 s by collecting samples at a resolution of 2-10 s per sample. 



Figure 2: Adsorption of reactants and factor Xa generation under flow. The reaction 
chamber was maintained at 37 C and perfused at 13.4 /i£/s with HEPES buffered medium, pH 
7.2, containing 0.15 N NaCl, 3 mM CaC£2, and 0.1 mM nonlabelled ovabumin. Maximal initial 
concentrations of reactant were 8 nM 3 H-fVIIa and 130 nM fXa. The TF density on the monolayer 
surface was estimated at > 25 fmol/cm 2 from both functional and immunological assays. (A) Total 
amounts of reactant (either fVIIa or fX, □) adsorbed to the monolayer were determined from the 
difference between the normalized concentrations of control tracer, 14 C(0)j an d test tracer, 3 H(a), 
collected in effluent samples. Tracer amounts were normalized as the fraction of the total added 
to the reaction chamber. (B) The amount of fX on the membrane (O) was determined from the 
difference between fX adsorbed (□) and fXa (•) released. The average rate of fXa generation for 
this experiment calculated from the slope of the middle linear segment of the progression curve 
(100-150 sec) was 2.1 ±0.02 fmol fXa/s/cm 2 . The mean from 13 similar experiments was 3.05 ±0.72 
fmol/s/cm 2 . 
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Figure 3: The lagged normal density curve and distributions of substrate and product 
concentrations. The open circles correspond to concentrations of fX (in nM) in aqueous phase 
determined from the concentration of nonreacting, nonadsorbing control tracer. The zero of the 
time axis is chosen to correspond to initial detection of 14 C. The qualitative fit to the lagged normal 
density curve (LNDC) yields the parameters <Ji ~ 21 s, r 2 ~ 91 s, and T2 ~ 19 s. The filled circles 
correspond to the fXa concentration released into each aliquot in the reacting system (xlO in the 
figure to facilitate comparison with fX values). 



Figure 4: Kinetic modeling of surface reactions. The set of nonlinear equations (Eqs. 0) 
was solved numerically and the associated parameters adjusted to obtain the best visual fit be- 
tween computed and experimentally measured fXa concentration curves. (A) Bulk fX concentration 
(O.lC^t), black) and product fXa (Tq, red). The approximate parameters achieving the best fit 
are: k+ = 15,k +a = 0.06, k_ a = Q,k +E = 0.0Q,k_ E = 0.0005, a lj2 = 0.8, a 6 = 0.1, fa = fa = 0.001, 
and (3q = 0.12. The amount of TF present in accessible enzyme complexes was assumed to be 0.32 
fmol/cm 2 . (B) Corresponding surface concentrations r2/100, T3,r4, and in fmol/cm 2 . 



Figure 5: Model predictions for alternative reaction mechanisms. (A) The predicted fXa 
generation as a function of time if the surface diffusion is much faster (k +a = 1.0) than that assumed 
in the simulations depicted in Fig. |], with all other parameters identical to those used in Fig. [|. 
(B) The corresponding surface concentrations. Note that membrane fX, T2, is much smaller than 
that in Fig| and estimated from measurements (not shown). (C) The predicted fXa production if 
the surface diffusion were slower than optimal. With k +a = 0.01, factor Xa is generated in lower 
quantities and at later times. (D) The predicted T 2 however, is much greater than that observed. 



Figure 6: Cartoon representing protein distribution on the cell surface. If the formation 
of E-fX occurs upon nearly each encounter of E and fX on the membrane surface, only E molecules 
near the perimeters of the domains will participate in catalysis of fX— >fXa. 
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TABLE I 

Reaction Chamber Parameters 



Dimensions (length x diameter) 0.8 cmxl.4 cm 

Volume 1.23 cm 3 

Microcarriers/chamber 6.7 ±0.2 x 10 4 

Microcarrier radius 7.5 x 10 -3 

Microsphere area 7.0 x 10 _4 cm 3 
Microcarrier/chamber volume fraction ~ 13% 

"Total reactive surface St ~ 47 cm 2 

6 TF antigen 0.63 fmol//ig protein 

TF activity 2.3 fmol//ig protein 

a Total surface estimated from microcarrier counts in samples from 3 different chambers. 
^Tissue factor antigen was measured by ELISA in detergent-lysed monolayers. Activity equivalent 
was estimated from TF activity in suspensions of intact and lysed (by one freeze/thaw cycle) cell 
monolayers, relative to the activity of standard dilutions of human recombinant TF reconstituted 
in phospholipid vesicles (30:70, PS/PC). The TF activity in suspensions of intact cell monolayers 
was 68.5 ± 19% of that in suspensions of lysed monolayers. 

TABLE II 

fXa generation rate as function of f Vila 



"Enzyme density 


b Rate 


c Flux(fmol/s) 




(fVIIa fmol/cm 2 ) 


(fmol fXa cm _2 s _1 ) 


fVIIa 


fX 


5.0-25 


2.9 ±0.05 


180 


3241 


2.2-9.0 


3.7 ±0.04 


52 


779 


1.9-3.4 


3.9 ±0.08 


23 


1728 


0.7-1.8 


3.7 ±0.06 


11 


608 










2000 



a Membrane density and adsorption rate of fVIIa were derived from the flux and adsorption rate 
coefficients previously measured (P9"D- Values are the initial and final concentrations measured 
during the linear interval of the progression curve. Tissue factor density in these experiments was 
estimated at > 25 fmol/cm 2 from cell protein and specific activity assays. 

b Average rate was calculated from the concentration of fXa measured by chromogenic assay in 
samples collected during the linear segment of the reaction progression curve. 
c Average fluxes of fVIIa and fX in each experiment were determined from the concentration of 
14 C-control tracer in the effluent. The total amounts of fVIIa and fX initially added to the 1.23 ml 
reaction chamber were 48-500 ng and 2800-9000 ng, respectively. 
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TABLE III 

Mean transit times of controls and fXa 



Factors 


Flow rate fA/s 


Mean transit time Tu 


(seconds) 






control 


fXa °, 


\ change 




21.4 ± 2.3 


56 ± 5 


79 ± 5.0 


41 ± 10 


a X, Vila 


11.7 ± 1.0 


107 ± 12 


144 ± 23 


32 ±7.0 




6.6 ±0.7 


247 ± 63 


270 ± 68 


12 ±3.0 




21.0 


54 


54 





fe Xa 


10.0 


109 


107 


-1 




8.0 


106 


167 


1 



a Mean transit times, Tp, of fXa and 14 C-ovalbumin were calculated using Eq. [| Factors Vila 
and X were added with 14 C-tracer, and fXa was generated in the chamber. Values are from four 
experiments at each flow rate. 

b Preformed fXa was added with 14 C-tracer. Values are from one experiment at each flow rate. 

TABLE IV 

Effects of flow rate on reaction 



J y (/d/s) a d{T 6 {t)/r 2 (t))/dt{xW 3 s- 1 ) b Yield (%) Average rate(fmol cm" 2 s" 1 ) 



4.6 


0.5 ±0.02 


17.6 


2.75 


6.0 


1.0 ±0.05 


17.5 


3.21 


6.4 


2.0 ±0.09 


17.1 


2.29 


13.3 


1.0 ±0.04 


7.8 


2.29 


14.2 


1.0 ±0.05 


5.2 


2.34 


27.0 


1.0 ±0.15 


6.0 


2.85 


28.0 


1.0 ±0.04 


3.9 


2.91 



a The fraction of membrane-bound fX converted to fXa increased linearly with time giving the 
slopes (roughly 0.001 s -1 ) indicated. The density of fX on the membrane was calculated from the 
difference between fX adsorbed and fX released as fXa. The concentration of the fXa released was 
measured directly in the effluent samples by amidolytic assay. 

b The yield of fXa is expressed as the percentage of the total fX added to the reaction chamber. 
Average reaction rates were determined from the steepest linear segment of the progression curves 
as in Fig. ||B. 
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Fig. 3 
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Fig. 5 
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